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We propose a hydrodynamics theory of collective quantum Hall states, which describes incom- 
pressible liquids, hexatic liquid crystals, a bubble solid and a Wigner crystal states within a unified 
framework. The structure of the theory is uniquely determined by the space-time symmetry, and 
a symmetry with respect to static shear deformations. In agreement with recent experiments the 
theory predicts two gapped collective modes for incompressible liquids. We argue that the presence 
of the above two modes is a universal property of a magnetized two-dimensional collective liquid. 
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I. INTRODUCTION 

Hydrodynamics or, more generally, continuum me- 
chanics is a common tool to describe collective degrees 
of freedom in condensed matter systems. This approach 
becomes especially powerful in the linear regime as it pro- 
vides a unified view on collective dynamics in solids, liq- 
uid crystals, normal liquids, superfluids, etc. (see, for ex- 
ample, Ref 1). The hydrodynamics equations of motion 
follow the universal conservation laws, while a particular 
state of matter enters the theory via a linearized "equa- 
tion of state" that is parametrized by a set of visco-elastic 
constants. A general form of the "equation of state" (nor- 
mally this is a dependence of the stress tensor on basic 
variables), and a number of independent parameters are 
determined by the macroscopic symmetry. 

Hydrodynamics represents the response of a many- 
particle system to external fields in terms of a few collec- 
tive variables. This assumes that for whatever reason 
the underlying particles loose their individual proper- 
ties and get collectivized. Apparently the above assump- 
tion is perfectly valid for quantum collective many-body 
states, such as superfluids or strongly correlated states 
of an electron gas in the fractional quantum Hall (FQH) 
regime 2 ' 3 . The hydrodynamics of superfluids has indeed 
a long and reach history (a general review as well as re- 
cent applications to trapped atomic Bose condensates can 
be found in Refs. 4,5). Yet the power of hydrodynam- 
ics approach in application to FQH fluids has been not 
fully explored. A notable exception is a low energy (in 
fact, adiabatic) extreme, where the structure of hydro- 
dynamic theory is well understood. In this regime the 
hydrodynamics of FQH liquids universally reduces to a 
topological Chcrn- Simons theory 6 (see also Ref. 7), which 
reflects a topological order of incompressible FQH states. 
An effective Chern- Simons theory successfully describes 
gapless edge excitations, but, by its low-energy nature, 
fails to reproduce fundamental bulk modes with a gapped 
spectrum 8 . 

A principal possibility to describe gapped magneto- 
roton modes within a continuum mechanics formalism 
has been demonstrated recently 9 . The most important 



observation of this work is that the gap A in the spec- 
trum of collective excitations requires a highly unusual, 
resonant frequency dependence of the shear stress. An ef- 
fective dynamic shear modulus should diverge at ui — A, 
otherwise collective modes will be gapless. It has been 
shown that the required resonant structure of the shear 
modulus naturally appears in hydrodynamics with an ad- 
ditional tensor collective variable that describes nontriv- 
ial precession dynamics of a shear stress. Using heuristic 
physical arguments in Ref. 9 we constructed the simplest 
equation of motion for the new tensor variable, and thus 
derived a "minimal" hydrodynamic theory that repro- 
duces a qualitatively correct dispersion of a gapped col- 
lective mode. 

The present paper reports a further progress in that 
direction. We show that there is a fundamental local 
symmetry behind the presence of a tensor field in the hy- 
drodynamics of incompressible FQH liquids. That is an 
invariance with respect to arbitrary static shear deforma- 
tions, which is a universal property of a liquid state of 
matter. The tensor dynamic variable is, in fact, a gauge 
field that supports the above local symmetry. 

From experimental side an interest to bulk collective 
excitations of FQH liquids has been recently renewed 
by the first observation of the long wavelength disper- 
sion of gapped modes at the filling factor v = 1/3. 10 
The most intriguing result of these experiments is a two- 
branch structure of the gapped excitation spectrum at 
small wave vectors. 

In the present paper we propose a universal hydro- 
dynamical interpretation of these experimental findings. 
We derive a hydrodynamics of FQH liquids, which re- 
lies solely on general symmetry arguments. The only as- 
sumption, which underlies our theory, is that the system 
is a two-dimensional (2D) collective liquid with a proper 
space-time symmetry (i. e. rotationally and time-reversal 
invariant). Incompressibility of the liquid, nonvanish- 
ing dissipationless Hall conductivity, a gapped spectrum 
of collective modes, etc., follow automatically from that 
assumption. The theory predicts two gapped modes 
(one with upward, and another one with downward dis- 
persions), which is in excellent agreement with recent 
experiments 10 . The key idea in constructing the hydro- 
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dynamics of incompressible liquids is to preserve the fun- 
damental property of any liquid state - a local symmetry 
under static volume-preserving deformations. By a suc- 
cessive reduction of the above local invariance we deduce 
theories of a hexatic liquid crystal, and two different crys- 
talline phases. 

The paper is organized as follows. In Sec. II we iden- 
tify basic symmetries, and present a phenomenological 
derivation of the linearized hydrodynamics of incom- 
pressible QH liquids. We construct a Lagrangian of the 
theory, derive equations of motion for basic collective 
variables, and analyze the spectrum of collective exci- 
tations. In Sec. Ill the theory is generalized to describe 
liquid crystal, and crystal states with hexagonal symme- 
try. In Sec. IV we summarize our main results. 

II. HYDRODYNAMICS OF INCOMPRESSIBLE 
QUANTUM HALL LIQUIDS 

A. Lagrangian of a magnetized collective liquid 

Our goal is to describe linear dynamics of an inter- 
acting 2D many-particle system confined to the x, y- 
plane, and subjected to a strong transverse magnetic field 
B = e z B. The hydrodynamic theory, which we derive in 
this section, is based on the following two, quite general 
assumptions about the ground/equilibrium state of the 
system: (i) The system is in a collective state that is ro- 
tationally, time-reversal 11 , and Galilee invariant; (ii) The 
system is in a liquid state. 

A few comments on the above assumptions are in or- 
der. The term collective state assumes the absence of a 
continuum of single particle excitations. Hence the dy- 
namics can be represented by a finite number of collec- 
tive variables. One of those variables, namely the dis- 
placement vector u(r, t), is fixed as it should appear in 
any linearized continuum mechanics. The velocity field 
v(r, t), and the variation of the density <5n(r,i) are re- 
lated to the displacement u as follows, 

v = d t u, Sn = — nVu, 

where n = v/2irl 2 is the equilibrium density (y is the 
filling factor, and I is the magnetic length). In general 
the displacement does not exhaust all relevant dynamic 
variables. It is well known that an adequate hydrody- 
namic description of ordered states of matter requires an 
introduction of additional fields that are specific for ev- 
ery system (a magnetization in ferromagnets, a director 
vector in nematic liquid crystals, a superfluid velocity in 
supcrfluids, etc.). We will see that a somewhat similar 
(though not precisely the same) situation takes place for 
incompressible FQH liquids. 

The assumption of a liquid state, (ii), can be also for- 
malized in terms of the displacement vector. By defini- 
tion, a liquid is a state of matter that does not respond 
to an arbitrary static shear (i. e. volume-preserving) de- 
formation. Mathematically shear deformations are de- 



scribed a divergenceless displacement vector. Hence by 
displacing a liquid system with a static, purely trans- 
verse vector f (r) we should not produce any physical ef- 
fect. Formally this means that the equations of motion 
for any 2D liquid should be invariant under the following 
transformation 12 

u' i {r,t)=u i {r,t)+Sijd j ^{r), (1) 

where ip(r) is an arbitrary function of only spatial vari- 
ables. The invariance under a volume preserving dif- 
feomorphism of Eq. (1) has a transparent interpretation 
within a Lagrangian formulation of fluid mechanics. It 
reflects the symmetry with respect to relabeling of dif- 
ferent infinitesimal fluid elements 13 . It is worth noting 
that the above diffcomorphism symmetry can serve as a 
most general mathematical definition of a liquid state. 
In terms of visco-elastic coefficients the symmetry under 
the transformation of Eq. (1) implies a vanishing static 
shear modulus. 

An important requirement, which we impose on 
the theory, follows from the universality of the high- 
frequency response. In the high-frequency, long wave- 
length limit any interacting many-body system behaves 
as an elastic continuum (see, for example, Ref. 14, and 
a recent paper of Ref. 15 for a most general demonstra- 
tion of this fact). Formally this means that in the high- 
frequency regime the exact stress tensor takes the stan- 
dard elastic form 16 

Pij = —SijKVu - n 00 {d t Uj + djUi - 5yVu), 

where the high-frequency shear modulus yUoo, and the 
bulk modulus K are universal constants that are related 
to the ground state properties 15 . Therefore in the limit 
u> — > oo the theory should reduce to the following equa- 
tion for displacement vector u, which describes linear dy- 
namics of the classical elastic medium 16 

md?u+Bd t uxe z -— V(Vu)-^V 2 u+VC/ H = F, (2) 

n n 

where J7h, arid F are the Hartree potential, and the ex- 
ternal force respectively. The high-frequency equation of 
motion, Eq. (2), is valid for all rotationally and Galilee 
invariant systems in the presence of an external mag- 
netic field. The above requirement uniquely determines 
the high-frequency form of the Lagrangian that gener- 
ates the correct left hand side of the equation of motion, 
Eq. (2) 

L°° = L + Lf, (3) 
where L is a linearized Lagrangian of an ideal liquid 17 

L = ^nm(d t u) 2 - ^nBu(d t u x e z ) - ^(Vu) 2 , (4) 

and Lg° is the standard elastic shear Lagrangian 

£s° = ~Hoo{diUjf. (5) 
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For the further discussion it convenient to introduce a 
traceless shear strain tensor 



Eq. (1): 



Sij = diUj + djUi - SijVu, 



(6) 



and to represent the shear Lagrangian of Eq. (5) in the 
following "natural" form 



Lg — - — /iooTrS _ —fJ^ooSij Sji, 



(7) 



One can straightforwardly check that Eqs. (5) and (7) 
are identical up to an irrelevant total derivative. It is 
important to outline that the high-frequency form of the 
Lagrangian, Eq. (3), is not an additional assumption, 
but the exact property of any rotationally and Galilee 
invariant interacting many-body system. 

Let us analyze the symmetries of the Lagrangian L°° 
defined by Eqs. (3), (4), and (7). Obviously it is rota- 
tionally, time- reversal, and Galilee invariant. [It is worth 
mentioning that the requirement of Galilee invariance 
fixes the coefficients in the first two terms in Eq. (4) to 
nm and nB respectively] Hence the our first assumption, 
(i), is fulfilled at the high frequency level. However, the 
assumption (ii), which corresponds to the local invari- 
ance of Eq. (1), is violated by the presence of the shear 
term in L°° . Apparently the first term, Lq, in Eq. (3) is 
invariant under the transformation Eq. (1). The second, 
shear term, Lg°, fails to satisfy this property since the 
shear strain tensor, Eq. (6), acquires a correction 

S'(r.f) = S(r,t) + <yS(r) ( 8 ) 

5S ij (Y)=£ ik d k d j %l)(v) + £ jk d k di%l>(v), (9) 

which changes Lg°. A common prescription to restore 
any local symmetry is offered by the gauge theory - one 
has to introduce a gauge field fj(r,t) that compensates 
unwanted changes in the Lagrangian. Namely, we replace 
S 2 in Eq. (7) by the form (S — r/) 2 , and add a Lagrangian 
L v that describes dynamics the gauge field f), i. e. 



rliq 



- iMl Tr(S 



flf+L. 



v 



The local symmetry of Eq. (1) is guaranteed if the term 
L v in the Lagrangian is invariant under the transforma- 
tion 



f7'(r,t)=fj(r,t)+5S(r), 



(10) 



where SS(r) is defined after Eq. (9). Additional restric- 
tions to the form of L g lq are imposed by the space-time 
symmetry, and a fixed high-frequency asymptotic form 
of Eq. (7). Apparently, since the shear strain tensor S, 
Eq. (6), is symmetric and traceless, the gauge field rj 
also has to be a symmetric, traceless tensor. Following 
the outlined route we construct the following most gen- 
eral (quadratic) shear Lagrangian that is rotationally and 
time-reversal invariant, and enjoys the local symmetry of 



4 iq - - 4^iTr(S - 17) 



1 



1 [i x 



4 A 2 



Tr(d t »)) 2 

+ A 2 Trr)(9 t r) x e z ) + 2A'TrS(9 t »7 x e z ) ,(11) 



where the cross product of e z and a symmetric tensor is 
defined as follows 

{rj x e z)ij = {sikVkj + £jkVki)/2- 

The term in the square brackets in Eq. (11) contains 
all allowed scalar invariants up to the second order in 
space-time derivatives. It should be noted that the rota- 
tional and the time-reversal symmetries allow two more 
quadratic invariants: 



Trr) 2 , Tr^r))^). 



(12) 



However these terms are forbidden by the invariance un- 
der the local transformation Eq. (10). As a result only 
the terms with at least one time derivative are left in L v . 

The final Lagrangian of a collective quantum Hall liq- 
uid takes the form 



L liq = L Q + 4 iq + nudt& - nuV<p 



(13) 



where L and Lg lq are given by Eqs. (4) and (11) respec- 
tively, and the last two terms describe the interaction 
with external potentials, a(r,i) and ip(r 7 t). 

The Lagrangian Eq. (13) defines a theory of two cou- 
pled fields - the vector field u(r, t) and a symmetric trace- 
less tensor field f)(r, t). Both our basic assumptions, (i) 
and (ii), are satisfied by Eq. (13). This is the most gen- 
eral harmonic Lagrangian for a 2D collective liquid in 
the presence of a transverse magnetic field. The theory 
contains four still undefined phcnomcnological parame- 
ters, /ii, A 2 , A2 and A'. The requirement of the high- 
frequency elastic response provides one constraint, which 
allows us to express /ii in terms of fi^. Three remain- 
ing parameters can be related to observable quantities 
- the excitation gaps, and the corresponding oscillator 
strength (see Sec. IIB below). 



B. Equations of motion and response functions 

The Lagrangian Eq. (13) generates the following equa- 
tion of motion for the displacement vector 

mndtU + nBd t u x e z - KV(Vu) + Vtt = nF, (14) 

where F = d t a — V<p, and 7r is the traceless shear stress 
tensor that is defined as follows 

7 r = - / i 1 [S-J7-(A7A 2 )a t r7xeJ (15) 

Similarly, variation of the Lagrangian with respect to 17 
yields the equation of motion for the gauge field 

5 2 r) + Ajrj - A 2 d t rj xe 2 = A 2 S + A'd t S x e z . (16) 
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The tensor field r) enters the equation of motion for 
the displacement vector u only via the stress tensor it, 
Eq. (15). Therefore it is more convenient technically to 
consider it as an independent dynamic variable. Making 
use of Eqs. (15) and (16) we can eliminate the gauge field 
r) in favor of it, and thus derive the following equation of 
motion for the shear stress tensor 

d 2 ir + Ajiz - A 2 d t ir x e z = -fn (l + A ,2 /A 2 ) <9 2 S 

-Mi(A 2 + 2A')d t S x e z . (17) 

Using Eq. (17) we can identify one of the phenomenolog- 
ical parameters entering the theory, namely the param- 
eter Hi. In the high-frequency regime the equation of 
motion for the displacement, Eq. (14), should reduce to 
its asymptotic form, Eq. (2). This means that at u> — > oo 
the solution to Eq. (17) should take the form it = — /XqoS. 
In the u — > oo regime the behavior of Eq. (17) is domi- 
nated by terms with the highest order of time derivatives, 
i. e. by the first terms in the right and left hand sides 
of Eq. (17). Thus the requirement of the elastic high- 
frequency response fixes parameter \x\ in the Lagrangian 
to the value 



m = Moo(l + A /2 /A 2 )- 1 . 



(18) 



The system of Eqs. (14), (17) completely determine the 
dynamics of the system. Using Eqs. (14) and (17) we can 
calculate the response to any configuration of external 
fields, and find all collective modes of the system. Let 
us first analyze the equation of motion for it, Eq. (17). 
This equation determines a dynamic "equation of state" 
that relates the stress tensor 7r to the strain tensor S. 
By symmetry the Fourier component of the shear stress 
tensor -k should be representable in the form 9 



it = -fi(u)S + iLuA(u)S x e z 



(19) 



where is the dynamic shear modulus, and A(o») is 

a "magnetic" modulus that is responsible for a Lorentz 
shear stress. Substituting the representation of Eq. (19) 
into Eq. (17) we arrive at the following system of two 
equations for functions /Lt(w) and A(w) 



(uj 2 - Aj)fi + u? A 2 A = uj 2 h c 



(cu 2 - A 2 )A + A 2A i = -fi c 



A 2 + 2A' 
1 + A' 2 /A 2: 



(20) 
(21) 



where we have used the relation of [i\ to fi^, Eq. (18). 
Straightforward solution of Eqs. (20), (21) yields the dy- 
namic moduli fi{w) and A(w) for the theory defined by 
the Lagrangian Eq. (13) 



A(w)=/i, 



/- 



uj 2 - A 2 



/-A- 



/+A+ 



j 2 - A 2 _ 



) 2 - A\ 



(22) 
(23) 



The quantities A + and A_ are related to parameters A 2 
and A 2 in the shear Lagrangian of Eq. (11) as follows 



A+A_=A 2 , A+-A_ = A 2 



(24) 



The oscillator strengths f± are given by the expressions 



f± = 



A±(A tT A') 



(A+ + A_)(A+A_ + A' 2 ) : 



(25) 



and satisfy the sum rule /++/_ = 1. 

The frequency dependence of shear modulus 
Eq. (22), agrees with our general expectations. It van- 
ishes at lu — * 0, and approaches fi^ at u> — > oo. A res- 
onant structure of Eqs. (22) and (23) is responsible for 
a gapped spectrum of collective excitations 9 (we will see 
that A± are, in fact, the excitation gaps). 

To calculate the response functions we need to sub- 
stitute the stress tensor, Eq. (19), into the equation for 
u, Eq. (14), and solve it for given external potentials 
a and ip. Below we consider the most interesting limit 
of a strong magnetic field when the first (acceleration) 
term in Eqs. (4) and (14) becomes irrelevant. Physically 
this corresponds the situation with all particles occupying 
the lowest Landau level (LLL). Formally the intra-LLL 
dynamics is described by the Lagrangian Eq. (13) [or, 
equivalently, by the equations of motions Eqs. (14), (17)] 
in the limit of vanishing bare mass, m — > 0. Substitut- 
ing Eq. (19) into Eq. (14), and neglecting the accelera- 
tion term we get the equation that describes the LLL- 
projected dynamics of the displacement vector u(q, u) 



-iu(nB + Aq 2 )u xe z + /j,q 2 u + ATq(qu) = nF, (26) 

where \i and A are given by Eqs. (22) and (23), and the 
force in the right hand side is related to the external 
potentials as follows, F = -iwa — iqep. 

Let us first calculate the Hall conductivity <th{u, q). 
By definition, (t#(cj,<7) relates the x-component of the 
current to the y-component of an electric field, 

jx = (?H(u,q)Ey 

Noting that j = nv = — iuonu. and solving Eq. (26) with 
a = 0, and — iqtp = E we find 



VH(u,q) 



u 2 (nB + Aq 2 ) 



LU 2 (nB + Aq 2 ) 2 - fi{p + K)q 4 ■ 



(27) 



The dissipationless (equilibrium) Hall conductivity is cal- 
culated from Eq. (27) by taking first the limit u> — * 0, and 
then the limit q — > 0. Since Eq. (22), vanishes as 

•J 2 , this order of limits gives the famous nonzero result 



(28) 



It is worth mentioning that in a normal viscous liquid as 
well as in a Fermi liquid at T = we have n(u>) ~ iui, 
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which leads to a vanishing dissipationless Hall conduc- 
tivity, a H (0,q) = 0. 

Using Eq. (26) one can easily show that the Hall con- 
ductivity also determines the response of the density, 
Sn = — nVu, to an external magnetic field b = V x a. In- 
deed, setting di(u>,q) = ieijqjb q /q 2 , and solving Eq. (26) 
we find 

Sn q = a H (io,q)b q 

Hence a nonvanishing limit in Eq. (28) implies a creation 
of a charge Q = v^/2tt by an adiabatic insertion of a 
magnetic flux <E> = b q= $. By inserting precisely one flux 
quantum, <f> = 27r, (the Laughlin's gcdanken experiment) 
we create a fractional charge v. 

Another important quantity is the density response 
function x(w,g), which relates Sn = — nVu to the ex- 
ternal scalar potential: Sn = x(w,g , )y>. The solution of 
the equation of motion, Eq. (26), with F = — iq<p yields 



n 



oj 2 (nB + Aq 2 ) 2 -fi(p + K)q 4 



(29) 



Inserting Eqs. (22) and (23) into Eq. (29) we find that 
at io —* and small wave vectors \ ~ <Z 4 > which signifies 
the proper incompressibility of the FQH liquid 8 . 

It should be stressed out that all abovementioned fun- 
damental low energy properties of FQH liquids are guar- 
anteed by the low frequency behavior of the shear mod- 
ulus, \i <~ lo 2 . In fact, this behavior guarantees that 
the low energy physics is dominated by the Lorentz force 
term in the Lagrangian/cquations of motion. In other 
words, in the limit lo — > the Lagrangian Eq. (13) re- 
duces to the form 



L liq ; 



nBu(d t u x e z ) + nud t a — nuS/tp 

nBsijUidtUj + nuidtcii — nuidiip, (30) 



which is precisely a Chern-Simons Lagrangian in a tem- 
poral gauge. Indeed, considering the standard Chern- 
Simons theory 



2,-1 I 

L = — — £ a p 1 A a d0A 1 - —e a p^A a dpa^, (31) 

47T I'K 

and setting Ai — 2imeijiij, Aq = 0, and ciq = lo we re- 
cover the low frequency Lagrangian Eq. (30). Thus our 
theory, defined by the Lagrangian of Eq. (13), smoothly 
interpolates between two exactly known limiting forms - 
the low energy Chern-Simons theory, Eq. (31), and the 
high-frequency elasticity theory, Eq. (3). Below we show 
that this interpolation predicts a spectrum of bulk collec- 
tive excitations, which is in excellent agreement with ex- 
perimental observations. It is worth noting that since in 
the low frequency limit our theory reduces to the Chern- 
Simons theory, it should reproduce the gapless low energy 
chiral edge excitations 6 . However, the full theory defined 
by Eq. (13) is also valid beyond the limit of asymptot- 
ically small frequencies. Therefore it should be able to 
predict a modification of the edge spectrum at higher 
energies. 



C. Spectrum of collective excitations 

To find collective modes of the system one has to set 
F = in Eq. (26), and to solve the resulting eigenvalue 
problem. Alternatively, the frequencies of eigenmodes 
can be determined from the poles of the density response 
function x(u,q), Eq. (29). The corresponding dispersion 
equation takes the form 

uo 2 [nB + A(uo)q 2 ] 2 - n(u) [ M (w) + K ] q A = 0. (32) 

Using the dynamic moduli /j,(u>), Eq. (22), and A(u>), 
Eq. (23), and solving Eq. (32) we find that the density 
response function has only two poles. These poles cor- 
respond to two collective modes with frequencies Q+(q) 
and Cl-(q). Both modes are gapped and have the follow- 
ing small-<? dispersion 



2,2 



(33) 



where /i^ = /loo/n is the high frequency shear modulus 
per particle, and I = l/\/B is the magnetic length. Since 
/± 5 Eq. (25), are positive, the dispersions of two collective 
modes always have opposite curvatures. 

The origin of gapped collective modes can be traced 
back to the dynamics of the tensor gauge field r), which 
is governed by the shear Lagrangian Eq. (11). In fact, 
the gaps, A + and A_, correspond to eigenvalues of the 
operator in the left hand side of Eq. (16). The presence of 
two modes with opposite curvatures of small-g dispersion 
relations agrees very well with recent experiments 10 . At 
q = both gapped modes are purely quadrupole, while at 
any finite q they have a mixed dipole-quadrupole char- 
acter with an increasing dipole component at larger q. 
Microscopically the two gaped modes are most likely re- 
lated to a mixture of small q magneto-roton excitations 8 
and two-roton bound states 18,19 . 

In the limit of small wave vectors the density response 
function, Eq. (29), takes the form 



x(io,q) = fi c 



f- 



+ 



L0 2 -n 2 _(q) L0 2 -Q 2 + (q) 



(34) 



Hence the quantities f± > determine the physical oscil- 
lator strengths for two collective modes. Thus Eqs. (24) 
and (25) relate all phenomcnological parameters of the 
theory to experimentally measurable quantities. 

Strictly speaking, the Lagrangian of Eq. (13) defines 
only a long wavelength theory. However, the results for 
collective modes look quite reasonable for any value of 
q. Indeed, a straightforward solution of Eq. (32) shows 
that with increase of q the function £1+ (q) monotonically 
increases, while Sl-(q) passes through a roton-like mini- 
mum, and approaches a constant value at ql — > oo. Thus 
an overall behavior of the lower branch, Q, -(q), surpris- 
ingly well reproduces a physically expected structure of 
the magneto-exciton spectrum, at least for the Laughlin 
fractions, v = I /(2k +1). 
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The most important result of the present theory is a 
prediction of a double-mode structure of gapped collec- 
tive excitations. The existence of two gapped modes is 
a generic fact for almost any set of parameters, A 2 , A 2 , 
and A' (it is somewhat more convenient to fix A + , A_, 
and A'). There is, however, a special, very narrow region 
in the parameter space where the upper branch of collec- 
tive excitations disappears, and only one gapped mode 
with a roton minimum is left. In this region our general 
construction reduces to a simplified single-mode magne- 
toelasticity theory of Ref. 9. In the single- mode theory 
the dynamic moduli, [i and A take the form 9 



uj 2 - A 2 ' 



A(w) = /i c 



' u 2 - A 2 ' 



(35) 



Using Eqs. (22)-(25) we find that expressions of Eq. (35) 
are recovered in two limiting cases: (1) A + — > oo, A_ = 
A, and A' is an arbitrary finite constant; (2) A + = A_ = 
A' = A. In this single-mode region of the parameter 
space an effective shear Lagrangian, Eq. (11), reduces to 
the following simple form 



riiq _ 



1 



= --MooTr(S 



<?) a + ^(^xe,). 



All the rest of the three-dimensional parameter space, 
(A + , A_, A'), corresponds to a theory with two gapped 
collective modes. It has been shown in Ref. 9 that the 
single-mode magnetoelasticity can be derived from the 
fermionic Chern-Simons theory 20 ' 21,22 ' 23 within the ran- 
dom phase approximation (RPA). In the present section 
we have demonstrated that the symmetry allows for a 
more general construction, which leads to a double-mode 
excitation spectrum. Microscopically this general FQH 
hydrodynamics should most likely emerge from post- RPA 
vertex corrections. 



III. HYDRODYNAMICS OF QUANTUM HALL 
CRYSTAL AND LIQUID CRYSTAL PHASES 

In Sec. II we derived a hydrodynamic theory of in- 
compressible quantum Hall liquids. In the present sec- 
tion this approach is generalized to non-liquid collective 
states. The hydrodynamics of FQH liquids was based on 
two principal assumptions, (i) and (ii). Note that at the 
long wavelength, i. e. for small q, in 2D systems the rota- 
tional symmetry is indistinguishable from the hexagonal 
one (with an accuracy up to q 2 ). Hence in the long wave- 
length limit the assumption (i) also covers all collective 
states of hexagonal symmetry. Let us keep the assump- 
tion (i), but relax (ii) by a successive reduction of the 
diffcomorphism symmetry defined by Eq. (1). 



A. Hexatic liquid crystal 

On the first stage we reduce the general local symmetry 
of Eq. (1) to the invariance under global shear deforma- 
tions. The later correspond to transformations Eq. (1) 



with <p(r) being a second order polynomial of r^. Such 
transformations produce a constant in space correction 
SS, Eq. (9), to the shear strain tensor. The rotational 
symmetry combined with the global invariance under the 
transformations Eq. (1) allow only one term to be added 
to the shear Lagrangian Eq. (11), 



6I#* = -jTr\;.,, . 



(36) 



where k is a positive constant. The full Lagrangian of 
a collective state that is invariant with respect to global 
shear deformations takes the form 



- hex 



L Q + 4 iq + <5£g ox + u<9 t a - uV<p. 



r hex 



(37) 



Apparently an additional term, <5Lg cx , in the La- 
grangian influences only the equation of motion for the 
gauge field f). Hence the equation for the displacement 
vector, Eq. (14), as well as the relation of Eq.(15) re- 
main unchanged. Only the equation of motion for the 
shear stress tensor -ft gets modified and acquires gradient 
corrections 



d^ + Aj 



1--V 2 
Mi 

+ A 2 kV 2 S - ^(A 2 + 2A')d t S x e 



7T - A 2 d t TT x e z = -fiood t S 

(38) 



[we note that Eq. (18) still holds, which follows from the 
universality of the u — > oo asymptotics] . 

The rotational, and the time-reversal symmetries 
require 9 the solution of Eq. (38) to be of the general 
form given by Eq. (19). Hence the strong field (m — > 0) 
expressions for the Hall conductivity and for the density 
response function, Eqs. (27) and (29), as well as the dis- 
persion equation, Eq. (32), do not change. The analytic 
structure of the elastic moduli \i and A is, however, differ- 
ent. Due to the gradient terms in Eq. (38) the dynamics 
moduli become functions of both the frequency ui and the 
wave vector q. 

The most important effect of the gradient term, <5Lg 0X , 
in the Lagrangian of Eq. (37) is a nonzero value of the 
shear modulus ji{uj, q) at u> = and finite q. Indeed, 
neglecting all time derivatives in Eq. (38) we get the fol- 
lowing static equation of state 



7T 



nq 



1 + k<7 2 //zi 



-m(o,<z)s. 



It the long wavelength limit the static shear modulus 
Ho{q) takes the form 



Ho(q) = n{0,q) w nq 2 



(39) 



(by construction the high-frequency limit of /i(w, q) is 
fixed to the value /ioo)- The static shear modulus Eq. (39) 
vanishes in the limit q — > 0. Therefore the state described 
by the Lagrangian Eq. (37) can not be a solid. To iden- 
tify the nature of this state we calculate the equal-time 
correlation function of the orientational order parameter 
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*(r) 



JVx 



u(r) 24 j^. g n j^ e temperature T the behav- 



ior of the correlation function (^*(r)W(0)) at r — > oo is 
determined by the static part of the Lagrangian L hcx . Us- 
ing an explicit form of Eq. (37) (without external fields), 
and performing the standard calculations we arrive at the 
following result 



(tf*(r)tf(0)) ~ 



(40) 



An algebraic decay of the correlator Eq. (40) is a clear 
signature of a hexatic liquid crystal 24 . 

A nonzero static shear modulus /io of hexatic state 
dramatically changes the low energy physics. Using the 
general formulas for the response functions, Eqs. (27) 
and (29), and the limiting form of the shear modulus, 
Eq. (39), we get 



lim lim <jH{u,q) = 0, 

q— >0 lo^O 



ft 

lim lim x(w,g) = - — . (41) 



Thus violation of the local symmetry of Eq. (1) destroys 
the dissipationless quantum Hall effect, and recovers the 
standard form of the compressibility sum rule 25 . 

The low energy behavior of the response functions, 
Eq. (41), is closely related to a dramatic modification 
of the structure of collective modes. In addition to two 
gapped modes, which remain practically unchanged in 
comparison with Eq. (33), we find one more nontrivial 
gaplcss solution to the dispersion equation, Eq. (32). The 
frequency of the gapless mode at small q is readily ob- 
tained by inserting the static shear modulus /zo, Eq. (39), 
into Eq.' (32): 



O = n Vmo(mo + K)q 2 l 2 

In a system with Coulomb interaction the small-g be- 
havior of the bulk modulus is dominated by the Hartree 
contribution, K w 2Trn 2 /q. 17 In this case the dispersion 
of the gapless mode in the hexatic phase takes the form 



(?) 



(for a short range interaction we get S!q cx <~ q 3 ). Inter- 
estingly, the above dispersion law is similar to that pre- 
dicted for a Goldstone mode in a quantum Hall nematic 
phase 26 . 

Recently a number of microscopic trial wave functions 
for different quantum Hall liquid crystal states have been 
proposed 27 . The single mode approximation (SMA) 8 ap- 
plied to the hexatic wave function yields a gapped col- 
lective mode with a roton minimum 28 . This behavior ex- 
cellently correlates with our results for the lower gapped 
mode, Cl-(q) - the dispersion of this mode has practi- 
cally the same form both for hexatic and for a liquid 
states. The microscopic SMA fails to produce a gapless 
mode that should be present in liquid crystal phases 28 . 
In contrast, our approach simultaneously predicts both 
the gapped modes and the dispersion of the Goldstone 
mode for the quantum Hall hexatic liquid crystal. 



B. Crystalline states 

Finally we completely destroy any kind of invariancc 
under shear deformations. The most radical way to do 
that is to drop out the gauge field r) by setting Ai, A 2 — > 
in Eq. (11). The resulting Lagrangian is simply L°°, 
Eq. (3). It describes the long wavelength dynamics of a 
hexagonal 2D solid with a frequency independent shear 
modulus /i = /loo. Apparently this state corresponds to 
a Wigner crystal. In the limit of a strong magnetic field 
(m — > 0) we get only one collective mode that is a well 
known 29 magneto-phonon with the dispersion 

= B- 1 ^IZq 3 / 2 . 

A more delicate way to violate the symmetry of Eq. (1) 
is to add to Eq. (37) the last invariant, Tr(j)) 2 , that is 
allowed by the rotational and the time-reversal symme- 
tries. The resulting total Lagrangian takes the form 



- bub 



- hex 



rM 2 Tr(j7) 2 , 



(42) 



where L hcx is defined after Eq. (37), and /J2 is an ad- 
ditional phenomenological constant. The second term 
in Eq. (42) leads to the following simple modification of 
the equation of motion for the shear stress tensor -k: in 
Eq. (38) the term — kV 2 is replaced by the combination 
-kV 2 + /i 2 . Noting this fact we easily realize that the 
Lagrangian of Eq. (42) describes a crystal state with a 
constant static shear modulus of the form 



fi Q = m(0,0) = 



M1M2 
Mi + M2 ' 



where \i\ is related to iioo by Eq. (18). 

A nonzero value of /x leads to a gapless magneto- 
phonon mode with a small-g dispersion similar to that 
for a classical Wigner crystal 



n 



^ h = B~ 1 ^2^ q 



3/2 



In addition, the dynamics of the gauge field 17, which is 
encoded in the equation of motion for the stress tensor 
■ft, still produce two FQH liquid-like gapped modes. The 
corresponding q = gaps, Ajj. ub and A?l ub , are related to 
parameters of the Lagrangian as follows 



A bub A bub = A 2 



1 + ^ 
Mi 



A bub _ A bub = ^ 



It is natural to interpret this intermediate state as a trian- 
gular lattice of liquid "bubbles" . The described structure 
of collective modes is indeed in a qualitative agreement 
with the Hartree- Fock excitation spectrum of a quantum 
Hall bubble phase 30 . 



IV. CONCLUSION 

In conclusion we proposed a unified hydrodynamical 
framework for the description of collective quantum Hall 
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states. In particular, the present theory covers incom- 
pressible FQH liquids, hcxatic liquid crystals, and two 
different hexagonal crystal phases that can be identified 
as a Wigner crystal, and a bubble solid. We predicted the 
dispersion of the Goldstonc mode for quantum Hall hex- 
atics, S1q cx <~ g 5 / 2 , and the presence of two gapped modes 
for incompressible quantum Hall liquids. The later re- 
sult naturally explains recent experimental observations 
of the long wavelength dispersion of collective modes in 
the v = 1/3 FQH state 10 . 

In short, the hydrodynamics of FQH liquids is based on 
the following observations. On the one hand, the high- 
frequency response of any system is universally elastic. 
On the other hand, a liquid state is insensitive to static 



volume preserving deformations. The former require- 
ment fixes the high-frequency form of the theory, while 
the later one implies a particular local symmetry which 
can be implemented using a proper gauge field. Since 
our approach relies only on symmetry, we believe that 
the predicted two gapped modes should be a universal 
property of a 2D collective magnetized liquid. 
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